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O 

Abstract 

a 

l — 1 We use multiple scattering and energy conservation arguments to 

describe dn/ drj\M A N A as a function of dn/drj\ pp in the framework of 
string percolation. We discuss the pseudo-rapidity r\ and beam rapidity 

f- — Y dependence of particle densities. We present our results for pp, Au- 

Au, and Pb-Pb collisions at RHIC and LHC. 

o 

\q 1 Introduction 

O 

As nuclei are made up of nucleons it is natural to look at nucleus-nucleus 
(A-A) collisions as resulting from the superposition of nucleon-nucleon (p-p) 
collisions, in the spirit of Glauber model approach and generalizations of it. 
X In the single scattering limit the average number of participating nucleons 

per nucleus, Na behave incoherently and 

dn dn 

^\n a n a = ^\ PP N A (1) 
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Eq (1) corresponds to the wounded nucleon model [lj[2j [3j. This model is 
expected to dominate at very low energy. In general, data do not agree with 
(1). 

At higher energy one has to take into account multiple scattering and 
one finds 

dn \ dn \ (aA+^s) at \ (o\ 

where 7Vi +a ^ is the estimated total number of nucleon-nucleon collisions 
and single scattering was subtracted [3]. 

It should be noticed that energy momentum conservation constrains the 
combinatorial factors of the Glauber calculus at low energy. The problem is 
that the energy momentum of Na valence strings has to be shared by N^ 3 
(mostly) sea strings. There are proposals to cure this problem, for instance, 
by reduction of the height of the rapidity plateau for sea strings [5]. In the 
same spirit, but reducing the effective number of sea strings rather than 
reducing the sea plateau, we write (see [1]). 

N f Ar(i+«M) (3) 

with 

a ^ = l^- TTi ( TTZd * (4) 

i> 1 + ln(y/s/so + 1) 

such that for yfs « -Jsq, a( v / s) — > 0, we are back to the wounded nucleon 
model, and for y/s » -y/so, a(y/s) — > |, and we have fully developed 
Glauber calculus. The need to take multiple scattering contribution was 
experimentally shown at RHIC [6j. 

Here as in [2] [3] our framework is the Dual Parton Model with parton 
saturation, and we work with Schwinger strings, with fusion and percolation 



In pp and Au-Au collisions or in general NaNa collisions the interactions 
occur with the formation of longitudinal strings in rapidity. The particle 
density dn / dy is expected to be proportional to the average number of strings 
(twice the number of elementary collisions) Nfa (see |3]), 



dn 
dy 



N A N A ~ N N A - ( 5 ) 



The string percolation model describes the multi particle production in 
terms of color strings stretched between the partons of the projectile and 
the target. In the impact parameter plane due to the confinement, the 
color of strings is confined to small area in transverse space Si = tttq with 
ro ~ .2 — .3 fm, these strings decay into new ones by qq — qq pair production 
and subsequently hadronize to produce the observed hadrons. In the impact 
parameter plane the strings appear as disks and as energy-density increases 
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the discs overlap, fuse and percolate, leading to the reduction of the overall 
color [8j[9j[10j. A cluster of n strings behaves as a single string with energy 
momentum corresponding to the sum of the individual ones. An essential 
quantity is the color reduction factor 



t \ 

N A < 



1 - e Vn a 

(6) 



where rf N is the string density in the impact parameter plane for NaNa 
collisions given as (see |U): 



V% a ^^N s Na , (7) 



Sn a is the area of the impact parameter projected overlap region of the inter- 
action covered by Na nucleons from nucleus A. Note that iVjy- = NpN^~ a 
and instead of (5) we have now 



dn 
dy 



NaN a ~ F^Na^Na ( 8 ) 



The color reduction factor F(rj N ) is a tool to slow down the increase of 
dn/dy with energy and number of participating nucleons. Note that in (1) 
nucleons interact incoherently and Sn a is in fact S p , while in (2), due to 
coherence, Sn a is the overall area of interaction. For details see We 
finally have at r] = 0, 



1 dn < dn \ h i F ^N A h A M^S) 1V1 , n s 



i 

N A dy^ A " A "dy™- ' F( V f 
with k being a normalization factor, 



and F(rf NA ) — > , t x and — > g, where iV® is the number of proton 

strings. At low energy N£ is around 2 growing with energy as e 2XY (faster 
than ^\pp) so that we can approximately write 

tf' = 2 + 4(£)V* y , (11) 
ftp 

We now generalize the results obtained in ref [3J. 

Based on the good description on data obtained by using the formula 
(9) for different atomic number and number of participants for different 
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energies at mid rapidity, we now apply the same formalism as used in pp 
to describe the rapidity evolution as suggested in ref [5j[llJ[12j obtaining a 
general formula for pseudo-rapidity dependence of A A collisions: 



1 dn ch 



N A drj 



/ i * 1 ' ' 

N A N A 

n exp ( s ) + 1 



i + A 



where J is the usual Jacobean J = coshri = anc [ K > _ K (exp ( M^ )- 

1). 

We now apply the formula to describe the charge multiplicity in p-p 
collisions for different energies in pseudo rapidity From our general formula 
(12) by using Na = 1 and A = 1, to consider p-p collisions the expression is 
reduced to 

K'nnbK—^^—- (13) 

n exp ( ( — s ) + 1 



drj 



2 Comparison with experimental data (RHIC, LHC) 

In figure (1) It is shown the comparison of the formula (13) applyied to 
different energies at different pseudo-rapidities with data from different ex- 
periments and energies, showing a good agreement in the evolution in pseudo 
rapidity and an increase in the plateau region as increasing with energy. 

In figure (2), (3), (4) and (5) it is shown the comparison between our 
results from formula (12) for Cu-Cu, Au-Au and Pb-Pb collisions at different 
energies, in agreement with data. 

In figure (6) We show some predictions for 3.2, 3.9 and 5.5 TeV energies 
at centrality — 5%, for Pb-Pb collisions. 

In the above computations we have used the following values of the pa- 
rameters: k = 0.63 ± 0.01, A = 0.201 ± 0.003, and = 245 ± 29 GeV, the 
same as obtained in [3], to describe the particle density ^|at a at a in the same 

power law as ^\pp- We had made here an extension to these descriptions to 
add the pseudo rapidity evolution with the same aim as in ref [5j. 

The new parameters values a ~ 0.34, 5 ~ 0.84, k\ = 1.2 had been set to 
adjust the equation (13) with data [13][14[15]. These results can be extended 
to describe proton- nucleus collisions. 
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Figure 1: Comparison of the results from the evolution of the dn c h/dr] with 
dependence in pseudorapidity from equation (13) for p-p collisions at differ- 
ent energies (lines), data is taken from ref. [13] |14| |15j . 

3 Conclusions 

We have discussed in a general way the physics of particle densities in pp and 
AA collisions. Our model gives a non-linear dependence of j^dn/dr)\AA on 
dn/dr)\pp. Particle densities, as a function of r\ and Y give a good description 
of Pb-Pb data at LHC in a wide region in r\. The same is observed for pp in 
a wide range of rapidity Y . 

Notice that recent data from TOTEM experiment measurements in the 
the charged particle pseudorapidity density dN c h/dr] in pp collisions at y/s = 
7 TeV for 5.3 < \tj\ < 6.4 have been compared to several MC generators and 
none of them has been found to fully describe the measurement, but our 
model is able to reproduce it. 
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Figure 2: Comparison of the results from the evolution of the dn c h/dr] with 
pseudorapidity from equation (12) for Cu-Cu collisions at 22.4 GeV, 62.4 
GeV, and 200 GeV energies, data is taken from ref. [16] .. Error bars in 
color blue, green, pink, red, purple and black are used for the corresponding 
centralities 45 - 55%, 35 - 45%, 258- 35%, 15 - 25%, 6 - 15%, - 6% 
respectively, lines in black show our results. 
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Figure 3: Comparison of the results from the evolution of the dn c h/dr] with 
pseudorapidity from equation (12) for Au-Au collisions at 19.6 GeV, 62.4 
GeV and 130 GeV energies, data is taken from [17] . Error bars in color blue, 
green, pink, red, purple and black are used for the corresponding centralities 
45 - 55%, 35 - 45%, 25 - 35%, 15 - 2§%, 6 - 15%, - 6% respectively, lines 
in black are the model results. 
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Figure 4: Comparison of the results from the evolution of the dn c h/dr] in 
pseudorapidity from equation (12) for Au-Au collisions at 200 GeV, data is 
taken from |17j . Error bars in color blue, green, pink, red, purple and black 
are used for the corresponding centralities 45 — 55%, 35 — 45%, 25 — 35%, 
15 — 25%, 6 — 15%, — 6% respectively, lines in black are the model results. 
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Figure 5: Comparison of the results from the evolution of the 
(dn c h/ drf) / {Np ar t/2) with the pseudorapidity from equation (12) for Pb-Pb 
collisions at 2.76 TeV, data is taken from |18j . Error bars in color green, blue, 
pink and red are used for the corresponding centralities 85 — 95%, 50 — 55%, 
— 90% and — 5% respectively, lines in black are the corresponding results 
from the model to the respective centrality, for the smaller centrality we use 
the number of participants corresponding to the 85-95% showed in dot and 
dashed line and in dashed line is the minimum number of participants equal 
2. 
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Figure 6: Predictions on the evolution of the {dn c h/ drj) / {N part /2) with pseu- 
dorapidity from equation (12) for Pb-Pb collisions at 3.2, 3.9 and 5.5 TeV 
energies at — 5% centrality. 
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